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Abstract 



^^^^ ^ We present simple analytical formulae which describe solutions of the RG equations 

h>. ' for Yukawa couplings in SUSY gauge theories with the accuracy of a few per cent. 

^O . Performing the Grassmannian expansion in these solutions, one finds those for all the 

^1 I soft couplings and masses. The solutions clearly exhibit the fixed point behaviour 

\^ ' which can be calculated analytically. A comparison with numerical solutions is made. 

o' 
a^ 

^ 1 Introduction 

Oh! 

r-| ! The renormalization group equations (RGEs) for the rigid couplings and soft parameters 

in SUSY gauge theories play a crucial role in applications. Actually, all predictions of the 
MSSM are based on solutions to these equations in leading and next-to-leading orders |l|. 
Typically, one has three gauge couplings, one or three Yukawa couplings (for the case of low 
or high tan/3, respectively) and a set of soft couplings. In leading order, solutions to the 
RGE for the gauge couplings are simple; however, already for Yukawa couplings, they are 
known in an analytical form only for the low tan/3 case, where only the top coupling is left. 
Moreover, even in this case solutions for the soft terms look rather cumbersome and difficult 
to explore 0. 

In a recent paper [^j , it has been shown that solutions to the RGE for the soft couplings 
follow from those for the rigid ones in a straightforward way.Q One takes the solution for the 
rigid coupling (gauge or Yukawa), substitute instead of the initial conditions their modified 
expressions 

(1) 
(2) 



ai =^ ai = 


= ai{l+M,r] + M,f] + 2MMvv), 


n ^ n = 


= Yk{l - AkV - AkV + AkAkVV + ^kVV) 
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where t] = 6'^,f] = 6"^, 9 and 9 are the Grassmannian parameters, and expand over these 
parameters. This gives the solution to the RGEs for the soft couphngs. 
Hereafter the following notation is used: 

«^ = t|L' y^^T^^ ^k = j:ml (3) 



167r^' IGtt^ 



=1 



where Qi and yk are the gauge and Yukawa couplings, respectively, and m'j are the soft masses 
associated with each scalar field. 

This procedure, however, assumes that one knows solutions to the ROE for the rigid 
couplings in the analytic form. For instance, in the case of the MSSM in the low tan/5 
regime this allows one to get solutions for the soft couplings and masses simpler than those 
known in the literature (see [0). At the same time, in many cases such solutions are 
unknown. Actual examples are the MSSM with high tan/3 and NMSSM. One is bound to 
solve the RGEs numerically when the number of coupled equations increases dramatically 
with the soft terms being included. 

Below we propose simple analytical formulae which give an approximate solution to the 
RGE for Yukawa couplings in an arbitrary SUSY theory with the accuracy of a few per cent. 
Performing the Grassmannian expansion in these approximate solutions one can get those 
for the soft couplings in a straightforward way. As an illustration we consider the MSSM in 
the high tan/3 regime. 

One can immediately see that approximate solutions obtained in this way possess infrared 
quasi-fixed points |^ which can be found analytically. They appear in the limit when the 
initial values of the Yukawa couplings are much larger than those for the gauge ones. Then, 
one can analytically trace how the initial conditions for the soft terms disappear from their 
solutions in the above mentioned limit. 

The paper is organized as follows. In Sect. 2, we consider the MSSM in the low tan/? 
regime, where all solutions are known analytically and describe briefly the Grassmannian 
expansion. In Sect. 3, we present our approximate solutions for the Yukawa couplings 
and obtain those for the soft terms. We also present numerical illustration and compare 
approximate solutions with the numerical ones. The fixed point behaviour is discussed. 
Section 4 contains our conclusions. The explicit formulae for the soft couplings and masses 
are given in Appendices. 

2 The MSSM: exact solutions in the low tan/? case 

Consider the MSSM in the low tan/3 regime. One has three gauge and one Yukawa coupling. 
The one-loop RG equations are 

33 

a- = -biaf, bi = {—,1,-3), i = 1,2,3. (4) 

5 

/'if} "1 Q \ 

Yt = Yt l^—a3 + 3a2 + —ai - 6Ytj . (5) 

with the initial conditions: aj(0) = ao, Yt{0) = Yq, and t = \n[Mx/Q'^). Their solutions are 
given by 



where 



Eit) = Y[il + kaotf^/^\ 



^15' ' 3^' 



(7) 



Fit) 



E{t')dt'. 



To get solutions for the soft terms, it is enough to perform the substitution a ^>- a and 
Y -^ Y for the initial conditions in (P) and expand over r/ and fj. Expanding the gauge 
coupling in (||) up to rj one has (hereafter we assume Mjo = '"^1/2) 



M,(t) 



"^1/2 



1 + hia^t 
Performing the same expansion for the Yukawa coupling and using the relations 
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mi/2itE-F). 
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one finds the well-known expression 
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To get the solution for the S term, one has to make expansion over 77 and 77. This can be 
done with the help of the following relations: 
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As a result one has 



E,(t)=^£zii + ^^°~ rny.QYoitE- F))' 



I+6Y0F 



QYoFY 



d^F 
drjdf] 


v,v 


2 ,2dE 


ml/2 


' d 
It 


ft'^dE\ 6F0 2dE' 
\E dt J 1+6^0^ dt 



(12) 



-JT • (13) 



which is much simpler than what one finds in the literature 0, though coinciding with it 
after some cumbersome algebra. 

One can also write down solutions for the individual masses using the Grassmannian 
expansion of those for the corresponding superfield propagators. For the first two generations 
one has 
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(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



where 
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The third generation masses get a contribution from the top Yukawa couphng 
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where A is related to S( (|T3D by 



Si — So 
So 



m 



1/2 



H2 



d /t^dE' 
~dt \E~dt 



m 



H, 



+ A/6, 
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(21) 
(22) 
(23) 
(24) 
(25) 
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With analytic solutions (^|TT],|T3]) one can analyze asymptotic and, in particular, find infrared 
quasi-fixed points which correspond to Yo ""^ oo 
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(26) 
(27) 
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One can clearly see that the dependence on Yq, Aq and So disappears from (p6[)-(p8D. Some 
residual dependence on ml is left for the soft masses and partially cancels with that of A. 

Below we demonstrate how the same procedure works in the case of approximate solu- 
tions. As a realistic example we take the MSSM in the high tan/3 regime. 



3 The MSSM: approximate solutions in high tan/3 case 

The one-loop RGE for the Yukawa couplings in this case look like 
Yt = Yt(^—as + 3a2 + —a,-6Yt-Y,y 

/'if} '~7 \ 

n = Yt(^—as + 3a2 + —a,-Yt-6Y,-Yrj, 
% = K (^3^2 + ^ai - 3n - 4k) . 

Since the exact solution is absent and might be too cumbersome, we look for an approximate 
one in a simple form similar to that of (^. 



3.1 The choice of approximate solution 



In choosing approximate solutions we follow the idea of f^ where an approximate solution for 
Yt and Yb ignoring 1^ has been proposed. Our suggestion is to consider separate brackets for 
each propagator entering into the Yukawa vertex. Then, one has the following expressions 
for the Yukawa couplings: 
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YboEb 
[1 + A{YtoFt + YboFb)]'^''[l + 2BYboFbf''[l + CiSY^oFb + Y^oFr)]'^"" 

YroEr 1 



^- 
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where the brackets correspond to the Q,U,H2, Q,D,Hi and L,E,Hi propagators, respec- 
tively. Here Ft and Ft are given by (|^) and (§) and Eb and Fr have the same form but with 



(r) 



cf> = (7/15,3, 16/3) and c\' = (9/5,3,0), respectively. 

The brackets are organized so that they reproduce the contributions of particular dia- 
grams to the corresponding anomalous dimensions. The coefficients A, B, C, A' and B' are 
arbitrary and their precise values are not so important. When Yukawa couplings Fjo ^-re 
small enough, one can make expansion in each bracket, and the dependence of these coeffi- 
cients disappears. However, for large couplings, which are of interest for us because of the 
fixed points, we have some residual dependence. The requirement that the sum of expo- 
nents equals 1 follows from a comparison with RGEs. Solutions are close to the exact ones 
when the brackets are roughly equal to each other. Apparently, since F^ < Ft ^ Fb and 
Yt- <Yb < Yt one cannot completely satisfy this requirement. Our choice of the coefficients 
A,B,C, A' and B' is dictated mainly by simplicity. In the following we choose them as 
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This gives approximate solutions like 
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(29) 
(30) 



(31) 



(32) 



(33) 



Solutions for Ai and Sj can be obtained by Grassmannian expansion with the initial condi- 
tions 

A(0) = v4o, S,(0) = So. (34) 

These initial conditions correspond to the so-called universality hypothesis which we follow in 
our numerical illustration for simplicity. However, one can choose arbitrary initial conditions 
for the soft terms when needed. It leads to an obvious modification of the formulae. 



One can get also the corresponding solutions for the individual soft masses. This can be 
achieved either by Grassmannian expansion of the corresponding brackets in (PT|)-(|5BD, or 
by expressing the masses through the Ss in an exact way. The second way gives a slightly 
better agreement with numerical solutions (see below). We present the explicit expressions 
for the soft terms and masses in Appendix A. 

3.2 Numerical analysis 

We start by investigating the precision of approximate solutions for the Yukawa couplings. 
To estimate the accuracy, we introduce a relative error which is defined as 

Y -Y 

^ approx ^numeric /or-\ 

^ = -^^ , (35) 

J- numeric 

and corresponds to the Mz scale {t = 66) at the end of the integration range. The accuracy 
for the solutions of soft terms is defined in the same way. 

Let us take at the beginning all three Yukawa couplings to be equal at the GUT scale 
and to have their common value Yq in the range (.01 -^ 25)ao- The upper limit is taken in 
order not to leave the perturbativity regime. We find that for Yq < ao the approximation 
errors are less than 3% for all Y's. While for Yt it remains smaller than 2% over the whole 
range of initial values at the GUT scale, for YJ, the error increases up to 4% and for 1^ up 
to 14% (for large values of Yq). It is worth mentioning that for small Yq (around ao/2 and 
below) the accuracy is very good (fractions of per cent or better). 

Consider now y^o = ^ro = ^o ^ lOao and let the top Yukawa coupling vary within 
the limits 1 < Yiq/Yq < 10 in order to examine the applicability of our formulae. In this 
case the accuracy it is spoilt a little bit with increasing initial values. Namely, the error 
for Yt increases up to 10%, and for Yj, and Yr up to 20%. However, if one keeps Yq in the 
range (1/100 -^ l/2)ao the accuracy for Yt remains better than 3%, and for Yj, and 1^ better 
than 10%. The particular case considered above seems to have the worst accuracy. This 
is not surprising since our approximate formulae are supposed to work best of all when all 
three Yukawa couplings are nearly equal. If we keep YrQ < Fw < Yto and the relative ratios 
less than 5, we get an average error of less than 5% for Yt, about 5% for Y^ and 10% for 
Yt-. This statement is illustrated in Fig.|I|. For each Yukawa coupling we have plotted the 
error as a function of YtQ in the range (1/10 4- 10)ao- The ratios are kept within the region 
1 < YtQ/YkQ < 5 and 1 < Fw/Ko < 3. 

Further on, we narrow the range of initial values up to (1/10 -^ 10)ao because the errors 
(defined as in (^5])) come to an asymptotic value for Yq > lOao and almost vanish for 
Yq < ao/10. The comparison of numerical and approximated solutions is shown in Fig.^ for 
three different sets of Yq^s. The approximate solutions follow the numerical ones quite well, 
preserving their shape, and they have a high accuracy, especially in the case of equal Yukawa 
couplings at the GUT scale. However, as can be seen from the top of Fig. |^, one can take 
arbitrary initial conditions for the Yukawa couplings, in particular those which are needed to 
fit the t/b/r masses, and to use our approximate solutions for these purposes. For the soft 
couplings. A, we take the initial values at the GUT scale to be Aq = (—2, —1, 0, 1, 2)mi/2 and 
leave Yqs in the narrow range as above. Then, we get an accuracy of (3 -i- 5)% for At and At,. 
For A-r the approximation is worse when Aq is taken to be negative or smaller than mi/2 (see 
Fig. H), but things go better for large initial values of Aq and we get an accuracy of about 
10%. Again it should be mentioned that this is an accuracy at the end point where A^- itself 
is close to and the accuracy defined as ( PB[ ) merely gives an odd hint of the precision. Along 
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Figure 1: Left side: Yf, Yj, and Y^- approximation errors. Numerical labels show the 
ratio Yho/YtQ for the corresponding groups of curves (split by different values of ^,-0). Right 
side: Comparison of approximate and numerical solutions for the soft masses. The curves 
correspond to the following choice of initial conditions: Aq = 0, ttIq = (l/2)m^ /2 and ml = 
2m\i2, YtQ = 5ao, Ybo = Fro = 2ao. Slepton masses are shown only for ml = m^^- Dotted 
lines correspond to analytical approximation, solid lines to numerical solution. 



the curves the accuracy is much better. In Fig.0 we have plotted the behaviour of At, Ai, and 
Ar for three different initial values of Aq, namely {— mi/2, 0, mi/2} and for one set of Yqs. As 
for the S's, keeping the range of parameter space for Yq and Aq as above, we get an accuracy 
of typically 2% for St (even better for fairly equal Yos). For E^ the precision is around 4%. 
With Hr we get into the same troubles as for Ar- The approximation becomes good (about 
10%) only for a large enough ratio of mQ/m^rt. The approximation errors for A's and S's 
are linked with those for Y. If one considers only the sets of small initial values for Yq (less 
than ao/2), then S's are approximated with a precision better than 1%, regardless of the Aq 
values. The precision for S increases with Aq, but this dependence is not so striking as the 
one on Yq. 

The approximate formulae for the soft masses may be derived from S using ( P^ ) 
In this case the approximate solutions give an accuracy of about 1^3% for m.Q^, m^ 
(3 -=- 5)% for m|), . For the Higgs masses we get a good approximation (of about 5% on 



Q3' "^^3 and 
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Figure 2: Comparison of approximate and numerical solutions for Y, A and S. Evolution 
of the soft terms A and S has been plotted for Fjo = 5ao ? ^w = ^to = 2ao- Dotted lines 
correspond to analytical approximation, solid lines to numerical solution. 



average) for m'jj^, and a satisfactory one for 'm'jj_^ (typically 10%). This accuracy is almost 
insensitive to the Aq variation (we took it to be in the range (—2 -^ 2)7711/2) and on the 
ratio rn1/mli2 (taken to be 0.5 -i- 2). The slepton masses (see Fig. |l]) are not approximated 
properly in an analogous way. This is mainly due to the less accurate approximation of Y^-. 
As a concluding remark on numerical analysis, it should be mentioned that one has a 
rather good approximation for small (less than ao/2) initial values of the Yukawa couplings. 
For larger values of F's one has a good approximation especially in the case of unification 
of the Yukawa couplings. 



3.3 The fixed points 



One can easily see that solutions (^)-(33) exhibit the quasi-fixed point behaviour when the 
initial values l^o = ^w = ^ro > Q^o- In this case, one can drop 1 in the denominator and the 
resulting expressions become independent of the initial conditions 
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FP 



Y, 



FP 



y 
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(36) 
(37) 
(38) 



These expressions being expanded over the Grassmannian variables give the quasi-fixed 
points for the soft terms and masses. The explicit expressions are presented in Appendix 
B. We see that the IRQFP behaviour is sharply expressed for Yt and 1^ (see Fig. |), and 
our approximate solution describes the fixed point line well. The same takes place for the 
corresponding As and Ss. For Y^^ A^ and S,- the accuracy is worse, however, the solution is 
still reliable. The soft mass terms there exhibit the same IRQFP behaviour, though some 
residual dependence on the initial conditions is left in full analogy with the exact solutions in 
the low tau/S case. The approximate solutions allow one to calculate the IRQFP analytically. 

One can see that the fixed points for the soft terms naturally follow from the Grass- 
mannian expansion of our approximate solutions (PB|)-(P5|) and they inherit their stability 
properties, as has been shown in |^. In particular, the behaviour of Ss essentially repeats 
that of the Yukawa couplings in agreement with |l^ . 

The existence of the IRQFPs allows one to make predictions for the soft masses without 
exact knowledge of the initial conditions. This property has been widely used (see, for 
example, |Tl[) and though the IRQFPs give a slightly larger top mass when imposing h — t 
unification, it is still possible to fit the quark masses within the error-bars and to make 
predictions for the Higgs and sparticle spectrum []12|. This explains general interest in the 
IRQFPs. 



4 Discussion 

We hope to convince the reader that the approximate solutions presented above reproduce 
the behaviour of the Yukawa couplings with good precision in the whole integration region 
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Figure 3: The IRQFP behaviour for Yukawa couphngs, soft terms Aj {, ,-,and S's. Numerical 
solutions for all A^s are done for Yto = 5ao, Vbo = ^ro = Sag. For S we took Aq = 0, three 
different sets of Yq's and ttIq. Dotted lines are the IRQFP from eqs. (p6| - p8D and Appendix B. 
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and for a large range of initial values. Relative accuracy is typically a few per cent and is 
worse only at the end of the integration region mainly due to the smallness of the quantities 
themselves. Moreover, we have shown how the approximate solutions for the soft terms and 
masses follow from those for the rigid couplings. This demonstrates how the Grassmannian 
expansion, advocated in , works in the case of approximate solutions as well. 

For illustration we have considered universal initial conditions for the soft terms. In recent 
time there appeared some interest in non-universal boundary conditions. Non-universality 



can also be included in our formulae at the expense of changing ([T0|) and (|3^) using the same 
substitution rules, see (|I]) and (§). 

Since the form of our approximate solutions has been "guessed" ad hoc starting from some 
reasonable arguments, there is no direct way to improve them. However, one can imagine 
more constructive derivation of those solutions which would allow one to make corrections. 
Needless to say that it is enough to construct a solution for the rigid terms. Solutions for 
the soft terms will follow automatically. 
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Appendix A 

We here present approximate expressions for the soft couplings and masses corresponding to 
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To find the individual soft masses one can formally perform integration of the RG equations 
and express the masses through Ss solving a system of linear algebraic equations. This gives 
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The masses of squarks and sleptons of the first two generations are given by (^)-(|T^. 



Appendix B 

We present here the IRQFPs for the soft couphngs and masses. They are obtained via 
Grassmannian expansion of (|36D -(|38D. 

.FP ( t dEt {tEt -Ft)+ [tEh -Fb) 5 {tEt - Ft)\ 



Ar^-^iM 



Et dt 7{Ft+Ft) 7 Ft J 

t dEb (tEt-Ft) + {tEh-Fb) 2{tEb-Fh) 3{tEh-Fb) + {tEr-Fr) 



Et dt 7{Ft + Ft) 7Fu 7{3Fk+Fr) 



pp ( t dE, A{tE,-Fr) 3{tE,-Fh) + {tEr-Fr 



^E^ dt 7Fr 7{3Fh+F^) 

s:FP ^ . f 2 t^'-^ + t^^ 5t^^ d (t^dEA 
* ^ '/'[ 7 (Ft + F,) 7 Ft ^dt[Et dt) 



2 [{tEt - Ft) + {tE, - n)]2 ^ 5 {tEt - Ft) 



7 (Ft + F,)^ 7 F^ ) 

Co J.2 dEt ,j.2dEj^ o j.2 dE^ o Of2 dE^ . ,2 dEr 
"^ ^ dt ^'^ dt "^ ^ dt '-' '-"^ dt ^'^ dt 

7 (Ft + F,) 7 F, 7 i3F, + F^) 

2 [{tEt - Ft) + jtE, - F,)]^ 2 {tE, - F,)' 
7 {Ft + F,)^ +7 F,2 

3 [3{tEb - Fb) + {tEr - F^)]^ d_ (^dEb\\ 
7 (3Fb + F,)2 ^ dt {Eb dt I) 

^ ~ ^i/2l^ 7 F, 7 {3F,+F,) ^dt[E,dt) 
^A{tEr-F,)^ , 3 [3{tE,-F,) + {tE,-F,)f\ 



7 F2 7 (3F, + F.)2 
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